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Introduction 

This work was inspired by the recent article 
"Mobius Transformations Revealed" by Dou- 
glas Arnold and Jonathan Rogness [3]. There 
the authors write: 

^^ Among the most insightful tools that mathe- 
matics has developed is the representation of a 
function of a real variable by its graph. . . . The 
situation is quite different for a function of a 
complex variable. The graph is then a surface in 
four dimensional space, and not so easily drawn. 
Many texts in complex analysis are without a 
single depiction of a function. Nor it is unusual 
for average students to complete a course in the 
subject with little idea of what even simple func- 
tions, say trigonometric functions, ^look like\^^ 

In the printed literature there are a few laud- 
able exceptions to this rule, such as the prize- 
winning "Visual Complex Analysis" by Tris- 
tan Needham [27j, Steven Krantz' textbook [20] 
with a chapter on computer packages for study- 
ing complex variables, and the Maple based 
(German) introduction to complex function the- 
ory [15] by Wilhelm Forst and Dieter Hoffmann. 
But looking behind the curtain, one encoun- 
ters a different situation which is evolving very 
quickly. Many of us have developed our own 
techniques for visualizing complex functions in 
teaching and research, and one can find many 



beautiful illustrations of complex functions on 
the internet. 

This paper is devoted to "phase plots", a spe- 
cial tool for visualizing and exploring analytic 
functions. Figure 1 shows such a fingerprint of 
a function in the complex unit disk. 




Figure 1: The phase plot of an analytic 
function in the unit disk 

The explanation of this illustration is deferred 
to a later section where it is investigated in de- 
tail. 

Phase plots have been invented independently 
by a number of people and it is impossible to 
give credit to someone for being the first. Origi- 



nally, they were mainly used in teaching as sim- 
ple and effective methods for visualizing com- 
plex functions. Over the years, and in partic- 
ular during the process of writing and rewrit- 
ing this manuscript, the topic developed its own 
dynamics and gradually these innocent illustra- 
tions transmuted to sharp tools for dissecting 
complex functions. 

So the main purpose of this paper is not only 
to present nice pictures which allow one to rec- 
ognize complex functions by their individual 
face, but also to develop the mathematical back- 
ground and demonstrate the utility and creative 
uses of phase plots. That they sometimes also 
facilitate a new view on known results and may 
open up new perspectives is illustrated by a uni- 
versality property of the Riemann Zeta function 
which, in the setting of phase plots, can be ex- 
plained to (almost) anyone. 

The final section is somewhat special. It re- 
sulted from a self-experiment carried out to 
demonstrate that phase plots are sources of in- 
spiration which can help to establish new re- 
sults. The main finding is that any meromor- 
phic function is associated with a dynamical 
system which generates a phase flow on its do- 
main and converts the phase plot into a phase 
diagram. These diagrams will be useful tools 
for exploring complex functions, especially for 
those who prefer thinking geometrically. 



Visualization of Functions 

The graph of a function f : D C C ^ C lives in 
four real dimensions, and since our imagination 
is trained in three dimensional space, most of 
us have diflSculties in "seeing" such an object j^ 

Some old books on complex function theory 
have nice illustrations of analytic functions. 
These figures show the analytic landscape of a 
function, which is the graph of its modulus. 




Figure 2: The analytic landscape of 

f{z) = {z-l)/{z' + z + l) 

The concept was not introduced by Johann 
Jensen in 1912 as sometimes claimed, but prob- 
ably earlier by Edmond Maillet [26j in 1903 (see 
also Otto Reimerdes' paper [31j of 1911). 
Differential geometric properties of analytic 
landscapes have been studied in quite a number 
of early papers (see Ernst Ullrich [34j and the 
references therein). Jensen [19j and others also 
considered the graph of |/p, which is a smooth 
surface. The second edition of "The Jahnke- 
Emde" JTE\ made analytic landscapes popular 
in applied mathematics. 

Analytic landscapes involve only one part of the 
function /, its modulus |/|; the argument arg/ 
is lost. In the era of black-and-white illustra- 
tions our predecessors sometimes compensated 
this shortcoming by complementing the analytic 
landscape with lines of constant argument. To- 
day we can achieve this much better using col- 
ors. Since coloring is an essential ingredient of 
phase plots we consider it in some detail. 
Recall that the argument argz of a complex 
number z is unique up to an additive multi- 
ple of 27r. In order to make the argument well- 
defined its values are often restricted to the in- 
terval (— 7r,7r], or, even worse, to [0, 27r). This 
ambiguity disappears if we replace argz with 
the phase zj\z\ of z. Though one usually does 
not distinguish between the notions of "argu- 
ment" and "phase" , it is essential here to keep 
these concepts apart. 

The phase lives on the complex unit circle T, 
and points on a circle can naturally be encoded 
by colors. We thus let color serve as the lacking 



^One exception is Thomas Banchoff who visualized four dimensional graphs of complex functions at [6]. 



fourth dimension when representing graphs of 
complex- valued functions. 





Figure 3: The color circle and the color-coded 
phase of points close to the origin 

The colored analytic landscape is the graph 
of I /I colored according to the phase of /. 
Since the modulus of analytic functions typi- 
cally varies over a wide range one better uses 
a logarithmic scaling of the vertical axis. This 
representation is also more natural since log |/| 
and arg/ are conjugate harmonic functions. 




Figure 4: The colored logarithmic analytic 
landscape of f{z) = {z — l)/(^^ + z + 1) 

Colored analytic landscapes came to life with 
easy access to computer graphics and by now 
quite a number of people have developed soft- 
ware for their visualization. Andrew Bennett 
[7j has an easy-to-use Java implementation, 
and an executable Windows program can be 
downloaded from Donald Marshall's web site 
[25] . We further refer to Chapter 12 of Steven 
Krantz' book [20], as well as to the web sites 
run by Hans Lundmark [21] and Tristan Need- 
ham |28j . Very beautiful pictures of (uncolored) 
analytic landscapes can be found on the "The 
Wolfram Special Function Site" [39jQ 



Though color printing is still expensive, colored 
analytic landscapes meanwhile also appear in 
the printed literature (see, for example, the out- 
standing mathematics textbook [T] by Arens et 
al. for engineering students). 
With colored analytic landscapes the problem 
of visualizing complex functions could be con- 
sidered solved. However, there is yet another 
approach which is not only simpler but even 
more general. 

Instead of drawing a graph, one can depict a 
function directly on its domain by color-coding 
its values, thus converting it to an image. Such 
color graphs of functions / live in the product 
of the domain of / with a color space. 
Coloring techniques for visualizing functions 
have been customary for many decades, for 
example in depicting altitudes on maps, but 
mostly they represent real valued functions us- 
ing a one dimensional color scheme. It is re- 
ported that two dimensional color schemes for 
visualizing complex valued functions have been 
in use for more than twenty years by now (Larry 
Crone [9], see Hans Lundmark [24j), but they 
became popular only with Frank Farris' review 
[T3] of Needham's book and its complement [14j . 
Farris also coined the name "domain coloring" . 
Domain coloring is a natural and universal sub- 
stitute for the graph of a function. Moreover it 
easily extends to functions on Riemann surfaces 
or on surfaces embedded M^ (see Konstantin 
Poelke and Konrad Polthier [30], for instance). 
It is worth mentioning that we human beings 
are somewhat limited with respect to the avail- 
able color spaces. Since our visual system has 
three different color receptors, we can only rec- 
ognize colors from a three-dimensional space. 
Mathematicians of the species gonodactylus oer- 
sted ijj could use domain coloring techniques to 
even visualize functions with values in a twelve- 
dimensional space (Welsch and Liebmann [37j 
p. 268, for details see Cronin and King [lOj). 
Indeed many people are not aware that natu- 
ral colors in fact provide us with an infinite 
dimensional space - at least theoretically. In 



^As of June 2009 Wolfram's tool visualizes the analytic landscape and the argument, but not phase. 
^ This is a species of shrimps which have 12 different photo receptors. 



reality "color" always needs a carrier. "Col- 
ored light" is an electromagnetic wave which is 
a mixture of monochromatic components with 
different wavelengths and intensities. A sim- 
ple prismatic piece of glass reveals how light is 
composed from its spectral components. Read- 
ers interested in further information are recom- 
mended to visit the fascinating internet site of 
Dieter Zawischa 




Figure 5: A typical spectrum of polar light 

The wavelengths of visible light fill an interval 
between 375 nm and 750 nm approximately, and 
hence color spectra form an infinite-dimensional 
space. 

How many color dimensions are distinguishable 
in reality depends on the resolution of the mea- 
suring device. A simple model of the human 
eye, which can be traced back to Thomas Young 
in 1802, assumes that our color recognition is 
based on three types of receptors which are sen- 
sitive to red, green, and blue, respectively. 
Since, according to this assumption, our visual 
color space has dimension three, different spec- 
tra of light induce the same visual impression. 
Interestingly, a mathematical theory of this ef- 
fect was developed as early as in 1853 by Her- 
mann Grassmann, the ingenious author of the 
"Ausdehnungslehre" , who found three funda- 
mental laws of this so-called metamerism [J£] 
(see Welsch and Liebmann [37j). 
Bearing in mind that the world of real colors is 
infinite dimensional, it becomes obvious that its 
compression to at most three dimensions cannot 
lead to completely satisfying results, which ex- 
plains the variety of color schemes in use for 
different purposes. The two most popular color 
systems in our computer dominated world are 
the RGB (GYM) and HSV schemes. 
In contrast to domain colorings which color- 
code the complete values f{z) by a two dimen- 
sional color scheme, phase plots display only 



f{z)/\f{z)\ and thus require just a one dimen- 
sional color space with a circular topology. As 
will be shown in the next section, they nev- 
ertheless contain almost all relevant informa- 
tion about the depicted analytic or meromor- 
phic function. 

In the figure below the Riemann sphere C (with 
the point at infinity on top) is colored using two 
typical schemes for phase plots (left) and do- 
main coloring (right), respectively. 





Figure 6: Golor schemes for phase plots and 
domain coloring on the Riemann sphere 

Somewhat surprisingly, the number of people 
using phase plots seems to be quite small. The 
web site of Francois Labelle [23j has a nice 
gallery of nontrivial pictureq^ including Euler's 
Gamma and Riemann's Zeta function. 
Since the phase of a function occupies only one 
dimension of the color space, there is plenty of 
room for depicting additional information. It 
is recommended to encode this information by 
a gray scale, since color (HUE) and brightness 
are visually orthogonal. Figure 7 shows two such 
color schemes on the Riemann t(;-sphere. 





Figure 7: Two color schemes involving 
sawtooth functions of gray 

The left scheme is a combination of phase plots 
and standard domain coloring. Here the value 
g of gray does not depend directly on log|i(;|. 



^Labelle justifies the sole use of phase by reasons of clarity and aesthetics. 



but is a sawtooth function thereof, g{w) = 
log \w\ — [log \w\\. This coloring works equally 
well, no matter in which range the values of the 
function are located. 

In the right scheme the gray value is the product 
of two sawtooth functions depending on log \w\ 
and w/\w\^ respectively. The discontinuities of 
this shading generate a logarithmically scaled 
polar grid. Pulling back the coloring from the 
t(;-sphere to the z-domain of / by the map- 
ping w = f{z) resembles a conformal grid map- 
ping^ another well-known technique for depict- 
ing complex functions (see Douglas Arnold [2j). 
Note that pulling back a grid instead of pushing 
it forward avoids multiple coverings. Of course 
all coloring schemes can also be applied to func- 
tions on Riemann surfaces. 
For comparison, the figure below shows the four 
representations of /(z) := {z—l)/{z'^ + z + l) in 
the square |Rez| < 2, |Imz| < 2 corresponding 
to the color schemes of Figure 6 and Figure 7, 
respectively. 





Figure 8: Four representations of the function 
f(z) = iz-l)/iz^ + z + l) 

Though these pictures (in particular the upper 
two) look quite similar, which makes it simple 
to use them in parallel, the philosophy and the 
mathematics behind them is quite different. We 
shall comment on this issue in the final section. 



The Phase Plot 

The phase of a complex function f : D ^ C 
is defined on Dq := {z e D : f{z) G C}, 
where C^ denotes the complex plane punctured 
at the origin. Nevertheless we shall speak of 
phase plots P : D ^ T,z ^ /(^)/l/(^)l on 
D^ considering those points where the phase is 
undefined as singularities. Recall that T stands 
for the (colored) unit circle. 
To begin with we remark that meromorphic 
functions are characterized almost uniquely by 
their phase plot. 

Theorem 1. // two non^zero meromorphic 
functions f and g on a connected domain D 
have the same phase, then f is a positive scalar 
multiple of g. 

Proof. Removing from D all zeros and poles of 
/ and g we get a connected domain Dq. Since, 
by assumption, f{z)/\f{z)\ = g{z)/\g{z)\ for aU 
z G Dq^ the function f/g is holomorphic and 
real- valued in Dq^ and so it must be a (posi- 
tive) constant. 

It is obvious that the result extends to the case 
where the phases of / and g coincide merely on 
an open subset of D. 

In order to check if two functions / and g with 
the same phase are equal, it suffices to compare 
their values at a single point which is neither a 
zero nor a pole. For purists there is also an in- 
trinsic test which works with phases alone: As- 
sume that the non-constant meromorphic func- 
tions / and g have the same phase plot. Then 
it follows from the open mapping principle that 
f ^ g \i and only if the phase plots of / + c 
and g^ c are different for one, and then for all, 
complex constants c 7^ 0. 



Zeros and Poles 

Since the phases of zero and infinity are unde- 
fined, zeros and poles of a function are singular- 
ities of its phase plot. What does the plot look 
like in a neighborhood of such points? 



If a meromorphic function / has a zero of degree 
n at Zq it can be represented as 

f{z) = {z-zorg{z), 

where g is meromorphic and g{zo) G C^. It fol- 
lows that the phase plot of / close to Zq looks 
like the phase plot of z^ at 0, rotated by the an- 
gle diTgg{zo). The same reasoning, with a nega- 
tive integer n, applies to poles. 




Figure 9: A function with a simple zero, a 
double zero, and a triple pole 

Note that the colors are arranged in opposite 
orders for zeros and poles. It is now clear that 
the phase plot does not only show the location 
of zeros and poles but also reveals their multi- 
plicity. 

A useful tool for locating zeros is the argument 
principle. In order to formulate it in the con- 
text of phase plots we translate the definition 
of winding number into the language of colors: 
Let 7 : T ^ L^o be a closed oriented path in 
the domain Dq of a phase plot P : Dq ^ T. 
Then the usual winding number of the mapping 
P o 7 : T ^ T is called the chromatic number 
of 7 with respect to the phase plot P and is 
denoted by chromp7 or simply by chrom7. 
Less formally, the chromatic number counts how 
many times the color of the point 7(t) moves 
around the complete color circle when 7(t) tra- 
verses 7 once in positive direction. 
Now the argument principle can be rephrased 
as follows: Let D be a Jordan domain with pos- 
itively oriented boundary dD and assume that 
f is meromorphic in a neighborhood of D. If f 
has n zeros and p poles in D (counted with mul- 
tiplicity)^ and none of them lies on dD, then 
n — p = chrom dD. 




Figure 10: This function has no poles. How 
many zeros are in the displayed rectangle? 

Looking at Figure 10 in search of zeros imme- 
diately brings forth new questions, for example: 
Where do the isochromatic lines end up? Can 
they connect two zeros? If so, do these lines 
have a special meaning? What about "basins of 
attraction"? Is there always a natural (cyclic) 
ordering of zeros? What can be said about the 
global structure of phase plots? We shall return 
to these issues later. 

The Logarithmic Derivative 

Along the isochromatic lines of a phase plot 
the argument of / is constant. The Cauchy- 
Riemann equations for any continuous branch 
of the logarithm log/ = In |/| + i arg/ imply 
that these lines are orthogonal to the level lines 
of 1/1, i.e. the isochromatic lines are parallel to 
the gradient of |/|. According to the chosen 
color scheme, we have red on the right and green 
on the left when walking on a yellow line in as- 
cending direction. 

To go a little beyond this qualitative result, we 
denote by s the unit vector parallel to the gradi- 
ent of I/I and set n := is. With Lp := arg/ and 
ip := log I/I the Cauchy-Riemann equations for 
log / imply that the directional derivatives of Lp 
and ip satisfy 

ds^ip = dn^ > 0, dn'ip = -ds(f = 0, 

at all points z of the phase plot where f{z) 7^ 
and f\z) y^ 0. Since the absolute value of dn^f 
measures the density of the isochromatic lines, 
we can visually estimate the growth of log|/| 
along these lines from their density. Because 



the phase plot dehvers no information on the 
absolute value, this does not say much about 
the growth of |/|. But taking into account the 
second Cauchy-Riemann equation and 

we obtain the correct interpretation of the den- 
sity dn^'- it is the modulus of the logarithmic 
derivative^ 



dn^ = 1/7/1 • 



(1) 



So, finally, we need not worry about branches 
of the logarithm. It is worth mentioning that 
dn^{z) behaves asymptotically like k/\z — Zq\ if 
z approaches a zero or pole of order fc at Zq. 
But this is not yet the end of the story. What 
about zeros of f'l Equation ([l]) indicates 
that something should be visible in the phase 
plot. Indeed, points Zq where /X^o) = and 
/(zo) 7^ are ^^color saddles'^ i.e. intersections 
of isochromatic lines. 

If /' has a zero of order /c at Zq, then z \-^ 
f{z) — f{zo) has a zero of order fc + 1 at Zq. 
Consequently / can be represented as 

f(z) = f{zo) + iz-zo)'^'giz) 

where g{zo) ^ 0. It follows that /(z) travels 
fc + 1 times around /(^o) when z moves once 
around Zq along a small circle. In conjunction 
with /(^o) 7^ this can be used to show that 
there are exactly 2/c + 2 isochromatic lines em- 
anating from Zq where the phase of / is equal 
to the phase of f{zo). Alternatively, one can 
also think of fc + 1 smooth isochromatic lines 
intersecting each other at Zq. 




Color saddles appear as diflFuse spots like in the 
left picture of Figure 11. To locate them pre- 
cisely it is helpful to modify the color scheme 
by superimposing a gray component which has 
a jump at some point t of the unit circle. If 
t := f{zo)/\f{zo)\ is chosen, then the phase plot 
shows a sharp saddle at the zero Zq of /' as in 
the right picture. 

Essential Singularities 

Have you ever seen an essential singularity? 
Here is the picture which usually illustrates this 
situation. 




Figure 12: The analytic landscape of 
fiz) = eV- 

Despite the massive tower this is not very 
impressive, and with regard to the Casorati- 
Weierstrass Theorem or the Great Picard The- 
orem one would expect something much wilder. 
Why does the analytic landscape not refiect this 
behavior? For the example the answer is easy: 
the function has a tame modulus, every contour 
line is a single circle through the origin. Now 
look at the phase plot in Figure 13: 




Figure 13: A phase plot depicting the essential 



Figure 11: Zeros of f are color saddles 



singularity of f{z) = e 



l/z 



But must there not be a symmetry between 
modulus and phase? In fact not. There is such 
a symmetry of modulus and argument (for non- 
vanishing functions) , but phase plots depict the 
phase and not the argument - and this makes a 
difference. 

So much for the example, but what about the 
general case? Perhaps there are also functions 
which conceal their essential singularities in the 
phase plot? 

In order to show that this cannot happen, we 
assume that / : D ^ C is analytic and has an 
essential singularity at Zq. 

By the Great Picard Theorem, there exists a 
color c G T such that any punctured neighbor- 
hood U of Zq contains infinitely many points 
Zk ^ U with f{zk) = c. Moreover, the set of 
zeros of /' in D is at most countable, and hence 
we can choose c such that c ^ /(^)/|/(^)| for 
all zeros z of /^ 

Since \c\ = 1, all isochromatic lines of the phase 
plot through the points Zk have the color c. As 
was shown in the preceding section, the modu- 
lus of / is strictly monotone along these lines. 
Now it is not hard to see that two distinct points 
Zk cannot lie on the same isochromatic line, be- 
cause these lines can meet each other only at 
a zero of /', which has been excluded by the 
special choice of c. 

Consequently any neighborhood of an essential 
singularity contains a countable set of pairwise 
disjoint isochromatic lines with color c. Com- 
bining this observation with the characteriza- 
tion of phase plots near poles and removable 
singularities we obtain the following result. 

Theorem 2. An isolated singularity Zq of an 
analytic function f is an essential singularity 
if and only if any neighborhood of Zq intersects 
infinitely many isochromatic lines of the phase 
plot with one and the same color. 

Note that a related result does not hold for 
the argument, since then, in general, the val- 
ues of SiTgf{zk) are different. For example, any 
two isochromatic lines of the function f{z) = 
exp(l/z) have a different argument. 



Periodic Functions 

Obviously, the phase of a periodic function is 
periodic, but what about the converse? 
Though there are only two classes (simply and 
doubly periodic) of nonconstant periodic mero- 
morphic functions on C, we can observe three 
different types of periodic phase plots. 




Figure 14: Phase plot of f{z) = e^ 




Figure 15: Phase plot of f{z) = sinz 




Figure 16: Phase plot of a Weierstrass 
p- function 

"Striped" phase plots like in Figure 14 always 
depict exponential functions f{z) = e^^+^ with 
a 7^ 0. Functions with simply p-periodic phase 
need not be periodic, but have the more gen- 
eral form e^^/^ g{z) with a G M and a p^periodic 



function g. Somewhat surprisingly, doubly peri- 
odic phase plots indeed always represent elliptic 
functions. 

The first result basically follows from the fact 
that the function arg / is harmonic and has par- 
allel straight contour lines, which implies that 
arg/(x + i^) = ax + /3y + ^. Since log |/| is 
conjugate harmonic to arg/, it necessarily has 
the form log |/(x + iy)\ = —ay + /3x + 5. 
If the phase of / is p-periodic, then we have 



h{z) := 



fiz+p) \fiz+p)\ 



m 



\f{z) 



and since h is meromorphic on C, it must be a 
positive constant e^. Now it follows easily that 
g{z) := f{z) • e~^^/^ is periodic with period p. 
Finally, if pi and p2 are periods of f/\f\ with 
Pi/P2 ^ K, then there exist ai, a2 ^ K such that 
f{z+pj) = e^^ /(^)- The meromorphic function 
g defined by g{z) := f'{z)/f{z) has only sim- 
ple poles and zeros. Integration oi g = (log f)' 
along a (straight) line from zq to Zq + Pj which 
contains no pole of g yields that 



an 



ZO+Pj 



g{z)dz. 



zo 



Evaluating now the area integral JJ^gdxdy 
over the parallelogram Q with vertices at 
0,Pi,P2,Pi + P2 by two different iterated inte- 
grals, we obtain a2Pi = aiP2- Since o^i, 0^2 ^ K 
and Pi/p2 ^ K this implies that aj = 0. 



Partial Sums of Power Series 

Figure 17 shows a strange image which, in sim- 
ilar form, occurred in an experiment. Since it 
looks so special, one could attribute it to a pro- 
gramming error. A moment's thought reveals 
what is going on here, at least at an intuitive 
level. This example demonstrates again that 
looking at phase plots can immediately provoke 
new questions. 




Figure 17: A Taylor polynomial of 
fiz) = 1/(1 -z) 

Indeed the figure illustrates a rigorous result 
(see Titchmarsh [33], Section 7.8) which was 
proven by Robert Jentzsch in 1914: 

// a power series Qq + aiz + a2Z^ + . . . has a pos- 
itive finite convergence radius R, then the zeros 
of its partial sums cluster at every point z with 
\z\ = R. 

The reader interested in the life and personality 
of Robert Jentzsch is referred to the recent pa- 
per [Hj by Peter Duren, Anne-Katrin Herbig, 
and Dmitry Khavinson. 

Boundary Value Problems 

Experimenting with phase plots raises a number 
of new questions. One such problem is to find 
a criterion for deciding which color images are 
analytic phase plots^ i.e., phase plots of analytic 
functions. 

Since phase plots are painted with the restricted 
palette of saturated colors from the color cir- 
cle, Leonardo's Mona Lisa will certainly never 
appear. But for analytic phase plots there are 
much stronger restrictions: By the uniqueness 
theorem for harmonic functions an arbitrarily 
small open piece determines the plot entirely. 
So let us pose the question a little differently: 
What are appropriate data which can be pre- 
scribed to construct an analytic phase plot, say, 
in a Jordan domain D7 Can we start, for in- 
stance, with given colors on the boundary dD7 



If so, can the boundary colors be prescribed ar- 
bitrarily or are they subject to constraints? 
In order to state these questions more precisely 
we introduce the concept of a colored set Kc^ 
which is a subset K of the complex plane to- 
gether with a mapping C : K ^ T. Any such 
mapping is referred to as a coloring of K. 
For simplicity we consider here only the fol- 
lowing setting of boundary value problems for 
phase plots with continuous colorings: 

Let D he a Jordan domain and let B he a con- 
tinuous coloring of its houndary dD. Find all 
continuous colorings C of D such that the re- 
striction of C to dD coincides with B and the 
restriction of C to D is the phase plot of an 
analytic function f in D. 

If such a coloring C exists, we say that the col- 
oring B admits a continuous analytic extension 
toD. 

The restriction to continuous colorings auto- 
matically excludes zeros of f in D. It does, 
however, not imply that / must extend contin- 
uously onto D - and in fact it is essential not 
to require the continuity of / on D in order to 
get a nice result. 

Theorem 3. Let D he a Jordan domain with 
a continuous coloring B of its houndary dD. 
Then B admits a continuous analytic extension 
to D if and only if the chromatic numher of B 
is zero. If such an extension exists, then it is 
unique. 

Proof. If (7 : D ^ T is a continuous coloring, 
then a simple homotopy argument (contract dD 
inside L^ to a point) shows that the chromatic 
number of its restriction to dD must vanish. 
Conversely, any continuous coloring B of dD 
with chromatic number zero can be represented 
as 5 = e^"^ with a continuous function ^p : dD -^ 
M. This function admits a unique continuous 
harmonic extension ^ to D. If ^ denotes a har- 
monic conjugate of $, then / = e^"^"^ is analytic 
in D. Its phase (7 := e^^ is continuous on D and 
coincides with B on dD. 



boundary colorings. This result can be gener- 
alized to phase plots which are continuous on 
D with the exception of finitely many singular- 
ities of zero or pole type in D. Admitting now 
boundary colorings B with arbitrary color index 
we get the following result: 
For any finite collection of given zeros with or- 
ders ni, . . . , nj and poles of orders Pi, . . . , p^ the 
houndary value prohlem for meromorphic phase 
plots with prescrihed singularities has a unique 
solution if and only if the (continuous) hound- 
ary coloring B satisfies 

chrom B = ni + . . . + Uj — pi — . . . — pk- 



The Riemann Zeta Function 

After these preparations we are ready to pay 
a visit to "Zeta", the mother of all analytic 
functions. Here is a phase plot in the square 
-40 < Rez < 10, -2 < Imz < 48. 




Theorem |3] parametrizes analytic phase plots 
which extend continuously on D by their 



Figure 18: The Riemann Zeta function 

We see the pole at z = 1, the trivial zeros at 
the points —2,— 4,— 6,... and several zeros on 
the critical line Rez = 1/2. Also we observe 
that the isochromatic lines are quite regularly 
distributed in the left half plane. 
Saying that Zeta is the mother of all functions 
alludes to its universality. Our starting point 
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is the following strong version of Voronin's Uni- 
versality Theorem due to Bagchi ^ (see also 
Karatsuba and Voronin [21j, Steuding [32]): 

Let D be a Jordan domain such that D is con- 
tained in the strip 



R 



{zgC : l/2<Rez< 1}, 



and let f be any function which is analytic in 
D, continuous on D, and has no zeros in D. 
Then f can be uniformly approximated on D by 
vertical shifts of Zeta^ Ct{^) •= C(^ + i^) ^ith 

teR. 

Recall that a continuously colored Jordan curve 
Jc is a continuous mapping C : J ^ T from a 
simple closed curve J into the color circle T. 

A string S is an equiva- 
lence class of all such col- 
ored curves with respect 
to rigid motions of the 
plane. Like colored Jor- 
dan curves, strings fall 
into different classes ac- 
cording to their chro- 
matic number. Figure 19 
depicts a representative 
of a string with chro- 
matic number one. 

We say that a string S lives in a domain D if 
it can be represented by a colored Jordan curve 
Jc with J C D. A string can hide itself in a 
phase plot P : L^ ^ T, if, for every ^ > 0, it has 
a representative Jc such that J C D and 




Figure 19: A 

representative of a 

string 



max \C{z) 



P{z)\ <e. 



In less technical terms, a string can hide itself 
if it can move to a place where it is invisible 
since it blends in almost perfectly with the back- 
ground. 

In conjunction with Theorem [3] the following 
universality result for the phase plot of the Rie- 
mann Zeta function can easily be derived from 
Voronin's theorem. 

Theorem 4. Let S be a string which lives in 
the strip R. Then S can hide itself in the phase 



plot of the Riemann Zeta function on R if it has 
chromatic number zero. 

In view of the extreme richness of Jordan curves 
and colorings this result is a real miracle. The 
three pictures below show phase plots of Zeta 
in the critical strip. The regions with saturated 
colors belong to R. The rightmost figure depicts 
the domain considered on p. 342 of Conrey's pa- 
per ^. 





Figure 20: The Riemann Zeta function at 
Imz = 171, 8230 and 121415 

What about the converse of Theorem [H? If 
there existed strings with nonzero chromatic 
number which can hide themselves in the strip 
i?, their potential hiding-places must be Jor- 
dan curves with non-vanishing chromatic num- 
ber in the phase plot. By the argument prin- 
ciple, this would imply that Zeta has zeros in 
R. If we assume this, for a moment, then such 
strings indeed exist: They are perfectly hidden 
and wind themselves once around such a zero. 
So the converse of Theorem [4] holds if and only 
if R contains no zeros of Zeta, which is known 
to be equivalent to the Riemann hypothesis (see 
Conrey [8j, Edwards [12j). 

Phase Flow and Diagrams 

Mathematical creativity is based on the inter- 
play of problem posing and problem solving, 
and it is my belief that the former is even more 
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important than the latter: often the key to solv- 
ing a problem lies in asking the right questions. 
Illustrations have a high density of informa- 
tion and stimulate imagination. Looking at pic- 
tures helps in getting an intuitive understanding 
of mathematical objects and finding interesting 
questions, which then can be investigated using 
rigorous mathematical techniques. 
This section intends to demonstrate how phase 
plots can produce novel ideas. The mate- 
rial presented here is the protocol of a self- 
experiment which has been carried out by the 
author in order to check the creative potential 
of phase plots. 

Let us start by looking at Figure 1 again. It 
depicts the phase of a finite Blaschke product^ 
which is a function of the form 






ZkZ 



with |z/e| < 1 and |c| = L Blaschke products are 
fundamental building blocks of analytic func- 
tions in the unit disc and have the property 
\f{z)\ = 1 for all z G T. The function shown in 
Figure 1 has 81 zeros Zk in the unit disk. 
Looking at Figure 1 for a while leaves the im- 
pression of a cyclic ordering of the zeros. Let us 
test this with another example having only five 
zeros (Figure 21, left). 




Figure 21: The phase plot of a Blaschke 
product with five zeros 

The left picture seems to confirm the expecta- 
tion: if we focus attention to the yellow color, 
any of these lines connects a zero with a certain 
point on the boundary, thus inducing a cyclic 
ordering. 

However, looking only at one specific color is 
misleading. Choosing another, for instance 



blue, can result in a diflFerent ordering. So what 
is going on here? More precisely: What is the 
global structure of the phase plot of a Blaschke 
product? This could be a good question. 
An appropriate mathematical framework to de- 
velop this idea is the theory of dynamical sys- 
tems. We here only sketch the basic facts; for 
details see [36j. 

With any meromorphic function / in a domain 
D we associate the dynamical system 



9i^) ■■- 



M f'iz) 



(2) 



The function g on the right-hand side of ^ ex- 
tends from Dq to a smooth function on D. This 
system induces a flow $ on L^, which we desig- 
nate as the phase flow of /. 
The fixed points of ^ are the zeros of / (re- 
pelling), the poles of / (attracting), and the ze- 
ros of /' (saddles). The remaining orbits are 
the components of the isochromatic lines of the 
phase plot of / when the fixed points are re- 
moved. Thus the orbits of the phase fiow endow 
the phase plot with an additional structure and 
convert it into a phase diagram. 
Intuitively, the phase fiow $ transports a col- 
ored substance ( "phase" ) from the zeros to the 
poles and to the boundary of the domain. 




Figure 22: Phase transport to the boundary, 
zeros of / and /' with invariant manifolds 

The left part of Figure 23 illustrates how 
"phase" of measure 27r emerging from the zero 
in the highlighted domain is transported along 
the orbits of $ until it is finally deposited along 
(parts of) the boundary. 

In the general case, where f : D ^ C is mero- 
morphic on a domain D and G C i^ is a Jordan 
domain with boundary J in Dq, the phase fiow 
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of any zero (pole) of / in G generates a (signed) 
measure on J. The result is a quantitative ver- 
sion of the argument principle which tells us in 
which way the phase of the zeros (poles) is dis- 
tributed along J (see Figure 23, right). 




Figure 23: The phase flow and the argument 
principle 

The question about the structure of phase plots 
of Blaschke products can now be rephrased in 
the setting of dynamical systems: What are the 
basins of attraction of the zeros of / with respect 
to the (reversed) phase flow? 
The key for solving this problem is given by 
the invariant manifolds of the saddle points, 
i.e., the points Qj G D where f'{aj) = and 




Figure 24: Invariant manifolds of the saddles 
and basins of attraction of the zeros 

Removing all unstable manifolds of the points 
aj from D results in an open set 5, which is 
the union of connected components Bj. Any 
component Bj contains exactly one zero bj of 
/, where multiple zeros are counted only once. 
The intersection of every set Bj with T is not 
empty and consists of a flnite number of arcs 
Aji. The complete set of these arcs covers the 
unit circle and two arcs are either disjoint or 
their intersection is a singleton. These sep- 
arating points are the endpoints of unstable 



manifolds which originate from saddle points. 
For later use we renumber the arcs Aji as 
^1,^2,...,^^ in counter-clockwise direction. 
It is obvious that the number s of separating 
points cannot be less than the number of dis- 
tinct zeros of /. In order to get an upper 
bound of s we assume that / has m distinct 
zeros with multiplicities /3i, . . . ,;5^ and k sad- 
dle points where f has zeros of multiplicities 
ai, . . . , a/e, respectively. Then we have 

ai + . . . + ak = m-1, /3i + ... + /3^ = n-l. 

The flrst equation follows from the well-known 
fact that the derivative of a Blaschke product of 
order n has exactly n — 1 zeros in D, this time 
counting multiplicity. From any saddle point aj 
exactly aj + 1 rays emerge which belong to the 
unstable manifold of aj. Since any separating 
point must be the endpoint of one such line, 
the total number s of separating points cannot 
be greater than k + /3i + ... + /3k = k + m— 1. 
Thus we flnally get 

m<s<m + k — 1. 

Examples show that both estimates are sharp. 
It turns out that the global topological struc- 
ture of the phase plot is completely character- 
ized by the sequence S of integers, which asso- 
ciates with any of the arcs yli, . . . , A^ (in con- 
secutive order) the number of the correspond- 
ing zero. This sequence depends on the spe- 
ciflc numbering of the zeros and the arcs, but 
an appropriate normalization makes it unique. 
For example, the Blaschke product depicted 
in Figure 24 is represented by the sequence 
^ = (1,2,3,2,4,5,4,2). 

Let us now return to Figure 1 again. Picturing 
once more that "phase" is a substance emerging 
from sources at the zeros which can exit the do- 
main only at its boundary, is it then not quite 
natural that phase plots of Blaschke products 
must look like they do? 

And if you are asking yourself what "natural" 
means, then this is already another question. 
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Concluding Remarks 

Phase plots result from splitting the informa- 
tion about the function / into two parts (phase 
and modulus), and one may ask why we do not 
separate / into its real and imaginary part. One 
reason is that often zeros are of special interest; 
their presence can easily be detected and char- 
acterized using the phase, but there is no way 
to find these from the real or imaginary part 
alone. 

And what is the advantage of using //|/| in- 
stead of In I /I ? Of course, zeros and poles can 
be seen in the analytic landscape, but they are 
much better represented in the phase plot. In 
fact there is a subtle asymmetry between mod- 
ulus and argument (respectively, phase). For 
example. Theorem [3] has no counterpart for the 
modulus of a function Jj 

Since phase plots and standard domain coloring 
produce similar pictures, it is worth mentioning 
that they are based on different concepts and 
have a distinct mathematical background. 
Recall that standard domain coloring methods 
use the complete values of an analytic function, 
while phase plots depict only its phase. Taking 
into account that phase can be considered as a 
periodization of the argument, which is (locally) 
a harmonic function, reveals the philosophy be- 
hind phase plots: Analytic functions are con- 
sidered as harmonic functions, endowed with a 
set of singularities having a special structure. 
Algebraically, phase plots forget about the lin- 
ear structure of analytic functions, while their 
multiplicative structure is preserved. 
This approach has at least two advantages. The 
first one is almost trivial: phase has a small 
range, the unit circle, which allows one visual- 
izing all functions with one and the same color 
scheme. Moreover, a one-dimensional color 
space admits a better resolution of singulari- 
ties. Mathematically more important is the ex- 
istence of a simple parametrization of analytic 
and meromorphic phase plots by their bound- 
ary values and their singularities (Theorem Isl). 



There is no such result for domain colorings of 
analytic functions. 

The following potential fields of applications 
demonstrate that phase plots may be a useful 
tool for anyone working with complex^valued 
functions. 

1. A trivial but useful application is visual in- 
spection of functions. If, for example, it is not 
known which branch of a function is used in a 
certain software, a glance at the phase plot may 
help. In particular, if several functions are com- 
posed, software implementations with different 
branch cuts can lead to completely different re- 
sults. You may try this with the Mathemat- 
ICA functions Log (Gamma) and LogGamma. An- 
other useful exercise in teaching is to compare 
the phase plots of exp(logz) and log(expz). 

2. A promising field of application is visual 
analysis and synthesis of transfer functions in 
systems theory and filter design. Since here the 
modulus (gain) is often more important than 
phase, it is recommended to use the left color 
scheme of Figure 7. 

3. Further potential applications lie in the area 
of Laplace and complex Fourier transforms^ in 
particular to the method of steepest descent (or 
stationary phase). 

4. Phase plots also allow to guess the asymp- 
totic behavior of functions (compare, for exam- 
ple, the phase plots of expz and sinz), and to 
find functional relations. A truly challenging 
task is to rediscover the functional equation of 
the Riemann Zeta-function from phase plots of 
( and r. 

5. Complex dynamical systems^ in the sense of 
iterated functions, have been investigated by by 
Felix Huang [17] and Martin Pergler [29] using 
domain coloring methods. The problem of scal- 
ing the modulus disappears when using phase 
plots; see the pictures of Francois Labelle [23j 
and Donald MarshaU [25]. 

6. The utility of phase plots is not restricted 
to analytic functions. Figure 25 visualizes the 



^ There is such a result for outer functions^ but it is impossible to see if a function is outer using only the boundary 
values of its modulus. 
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function 

h{z) := Im (e-T z^) + i Im (e^ {z - 1)^) , 

with n = 4. This is Wilmshurst's example j38] 
of a harmonic polynomial of degree n having the 
maximal possible number of n^ zeros. For back- 
ground information we recommend the paper on 
gravitational lenses by Dmitry Khavinson and 
Genevra Neumann l22l. 




Figure 25: A modified phase plot of 
Wilmshurst's example for n = 4 

To understand the construction of the depicted 
function it is important to keep track of the ze- 
ros of its real and imaginary parts. In the figure 
these (straight) lines are visualized using a mod- 
ified color scheme which has jumps at the points 
1, i, — 1 and — i on the unit circle. 

Besides these and other concrete applications 
one important feature of phase plots is their 
potential to bring up interesting questions and 
produce novel ideas. If you would like to try 
out phase plots on your own problems, you may 
start with the following self-explaining Matlab 
codeH 

xmin=-0 . 5 ; xmax=0 . 5 ; ymin=-0 . 5 ; ymax=0 . 5 ; 

xres = 400; yres = 400; 

x = linspace(xmiii,xmax,xres) ; 

y = linspace(ymiii,ymax,yres) ; 

[x,y] = meshgrid(x,y) ; z = x+i*y; 

f = expCl./z) ; 



p = surf (real (z) ,imag(z) ,0*f, angle (-f)) ; 
set (p , ^ EdgeColor ^ , ^ none ^ ) ; 
caxisC [-pi,pi] ) , colormap hsv(600) 
view(0,90), axis equal, axis off 

Though the phase of a function is at least of the 
same importance as its modulus, it has not yet 
been studied to the same extent as the latter. 
It is my conviction that phase plots are problem 
factories^ which have the potential to change 
this situation. 

Technical Remark. All images of this article 
were created using Mathematica and Mat- 
lab. 
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